Abstract 



We are interested in formulating a viscous model of the universe based on The 
Bianchi Type IV algebra. We first begin by considering a congruence of fluid lines in 
spacetime, upon which, analyzing their propagation behaviour, we derive the famous 
Raychaudhuri equation, but, in the context of viscous fluids. We will then go through 
in great detail the topological and algebraic structure of a Bianchi Type IV algebra, 
by which we will derive the corresponding structure and constraint equations. From 
this, we will look at The Einstein field equations in the context of orthonormal frames, 
and derive the resulting dynamical equations: The Raychaudhuri Equation, generalized 
Friedmann equation, shear propagation equations, and a set of non-trivial constraint 
equations. We show that for cases in which the bulk viscous pressure is significantly 
larger than the shear viscosity, this cosmological model isotropizes asymptotically to 
the present-day universe. We finally conclude by discussing The Penrose-Hawking 
singularity theorem, and show that the viscous universe under consideration necessarily 
emerged from a past singularity point. 
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1 Introduction 



For centuries past and into the present time, much of physics has been taught and done 
in the flat, Euclidean space, with time being considered as a dependent, constant vari- 
able. Indeed, two such fields of physics have emerged. There is a local physics, which 
involves formulating theories of local phenomena such as fluid dynamics, quantum me- 
chanics, and electromagnetism. On the other hand, there is an universal physics which 
attempts to describe the nature of the universe on a large-structure scale. While the 
former have been described extremely successfully by The Navier-Stokes equations, 
Schrodinger's Equation, and Maxwell's equations, only Einstein's Theory of General 
Relativity has been able to successfully, and completely describe the latter. In recent 
times, many attempts have been made with respect to combining the two into one grand 
perspective of describing the whole universe, from the very small to the very large. Ein- 
stein's theory conveniently allows for this by the notion that matter induces curvature 
in the spacetime manifold, and the resulting spacetime geometry prescribes the mo- 
tion of of this matter. In the context of this description for how our universe really 
works, we attempt to formulate a theory of viscous fluids in a spatially homogeneous 
spacetime. This is not as easy as simply applying the Navier-Stokes equations. Despite 
their overwhelming success in the classical theory of fluid dynamics, the Navier-Stokes 
equations only apply in Euclidean space. The main question is that of the classifi- 
cation of partial differential equations. In theory, since the Navier-Stokes equations 
are a coupled set of hyperbolic partial differential equations, and therefore, constitute 
causal solutions, one can map these over to the pseudo-Riemannian manifold of Gen- 
eral Relativity by the usual partial to covariant derivative rule. However, in the case of 
viscosity, the Navier-Stokes equations become elliptical partial differential equations, 
and therefore lead to acausal solutions. This is not allowed in any relativistic theory, 
since this would technically imply that fluid propagation speeds are able to exceed the 
speed of light. Therefore, other methods must be used to describe the viscous fluid flow 
in a General Relativistic context. We are particularly interested in viscous/imperfect 
fluids as compared to inviscid/perfect fluids for the sole reason that inviscid flow is 
an idealization. In the context of cosmology, many astrophysical fluids such as gases 
exhibit high compressibility, where viscosity must be included in any fluid analysis for 
more realistic results. In addition, to the best of the author's knowledge at the time 
of the writing of this thesis, viscous flows in the Bianchi Type-IV universe have not 
been described, so this thesis also is an opportunity to present some new and valuable 
information. Perhaps the most important application is to early-universe cosmology. 
The early universe consisted of high-temperature fluids, which necessarily means that 
viscosity must be included to appropriately describe any fluid matter. In addition, the 
early universe is taken to be spatially homogeneous and anisotropic, which somehow 
isotropized to the current FLRW models of the present-day universe. 

We will proceed in the following manner: We will first derive the viscous analogue 
of the all- important Raychaudhuri equation which describes the expansion/contraction 
behaviour of a congruence of fluid lines, and necessarily the universe. After discussing 
some qualitative solutions to this equation, we solve The Einstein Field equations for 
a Bianchi Type-IV algebra with a viscous fluid matter source using the approach of 
orthonormal frames. The point of this approach is that the field equations which typi- 
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cally constitute ten coupled, non-linear, hyperbolic partial differential equations reduce 
to a system of autonomous, first-order, coupled ordinary differential equations which 
are easier to solve using numerical methods. We will conclude with a discussion of 
The Penrose-Hawking singularity theorem applied to a viscous universe. We assume 
that the reader has adequate knowledge of general relativity and differential geometry 
throughout. In addition, we use geometrized units, such that c = G = 1. In ad- 
dition, as typical in cosmology settings, we will use comoving coordinates, such that 
u a = (1, 0, 0, 0), the metric signature (—1, 1, 1, 1), which implies that u a u a = — 1. 

2 The Notion of Spatial Homogeneity 

Standard cosmological models have for the most part taken the universe to be at 
a minimum, spatially homogeneous. The historical reasons for this have been well- 
described in the literature. The interested reader should see for example [41] and 
[24]. In terms of the timeline of our universe, cosmological models are essentially 
divided up into three categories. The first category consists of cosmological models that 
describe the immediate post-Big Bang universe, that is the universe on Planck scales. 
Such models are the typical quantum gravity models based on the concept of Misner's 
minisuperspace [26] [28] and The Hamiltonian formulation of General Relativity given 
by Arnowitt, Deser, and Misner [1] , and Dirac [10]. There are of course others that 
have contributed to this subject, however, the main work seems to be accredited to 
ADM and Dirac. The second category consists of models of the early universe, that 
is, the post-Big Bang universe considered on scales much larger than Planck scales 
such that classical General Relativity can be applied. The key point of these models is 
that they are taken to be spatially homogeneous, but since it is widely assumed that 
this stage of the universe still had much anisotropy in terms of cosmic radiation, these 
models are taken to be anisotropic, and hence conveniently characterized by the Bianchi 
classifications. Much work has been done in this area as well. For example, Singh and 
Baghel have recently studied Bianchi Type V bulk viscous cosmological models with 
a time dependent cosmological constant. [36] Misner famously studied the Bianchi 
Type IX/Mixmaster models which describe a chaotic decay of the anisotropy of the 
universe. [25] In addition, Belinskii, Khalatnikov, and Lifshitz have studied in detail 
the oscillatory approach to a singular point in the relativistic cosmology. [3] The third 
category consists of the present-day models of the universe which are entirely described 
by The Friedmann-LeMatire-Robertson- Walker (FLRW) models with a perfect fluid 
matter source. Not only are these models spatially homogeneous, but they are also 
isotropic about every point making them a special class of the more general Bianchi 
models. One of the major areas of research in cosmology is to understand how the 
anisotropies of the early universe decayed to the present-day state of isotropy. This 
problem has been looked at in quite some detail as well. The interested reader is 
encouraged to read [16], [4], [2], especially [20], [9], and [27]. In addition, Salucci and 
Fabbri [34] have in great detail analyzed the cosmological evolution of general Bianchi 
models in the adiabatic regime. They showed that a constant or slowly varying energy- 
momentum tensor isotropized the Hubble expansion. As alluded to earlier, our area 
of interest for this thesis is to develop a model of the early universe with a viscous 
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fluid matter source. This model is unique in that the Bianchi Type-IV model we will 
consider has not been described in the literature. We have so far said a great deal 
about spatial homogeneity and (an)isotropy. Therefore, it is appropriate to describe 
what these terms mean in a more precise, mathematical context. 

2.1 The Mathematical Foundations of Spatial Homogeneity 

In this section, we will build the mathematical foundations for spatial homogeneity. 
The arguments we present here are largely based on the descriptions given by Hervik 
and Gr0n [15], Ellis and MacCallum [14], Eisenhart [12] [13], and Stephani et.al [37]. 
The interested reader is encouraged to refer to these sources for further details. How- 
ever, we attempt to provide as comprehensive of a description as possible. Basically, a 
space is homogeneous if it admits a group of motions. That is, a space is topologically 
homogeneous if you can carry one point to any other point via an isometry. If we 
consider a space M on the manifold described by the metric tensor g, we can then 
define an isometry group as: 

Isom(M) = {(f) : M M\tj> isometry} (2.1) 

An isometry operation on the metric is one that leaves the metric unchanged, that 
is, <f> * g = g- From classical Riemannian geometry, we know that Killing vector fields 
generate isometries. We can therefore say that a space is homogeneous if for each 
pair of points (p,q) & M, 3(j> € Isom(M) such that (j>(p) = q. In some literature [22], 
homogeneous spaces are referred to as transitive. 

We can extend this definition a bit further. We could define an isotropy subgroup 
of a point p € M by: 

S P (M) = {<j> e Isom(M)\<t>(p) = p} (2.2) 

Intuitively, this definition is providing for a mechanism for spatially homogeneous 
spacetimes that are also isotropic. By leaving one point fixed as the above definition 
implies, we are allowing for the possibility of rotations. We shall see the importance 
of this subgroup in due course. 

If we further let dim{I som(M)} = n, and dim{S p (M)} = m, a space is homo- 
geneous if n > dim{M}. It is said that M is simply transitive if it is homogeneous 
and n — dim{M}. M is said to be multiply transitive if M is homogeneous, and 
n > dim{M}. 

In addition, if we consider for a point p € M in a subspace of M such that: 

H p = {q£ M\q = <f>(p) for a G Isom(Al)} (2.3) 

This subspace is known as the orbit of p under the isometry group. Going back 
to the idea of moving from one point to another in a space, all points that we can 
reach via the action of the isometry group on p, is necessarily in the orbit of p. In 
addition, if the orbit of p is actually the whole space, that is, H p = M, then space is 
homogeneous/transitive. 

In a simply transitive space, there always exists a basis of Killing vectors such 
that: 
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= (2-4) 

This indeed is the main method described by Stcphani [37] and Taub [38]. However, 
this seems slightly inefficient as to get an expression for the metric tensor involves 
expressing these Killing vector basis vectors in canonical form, that is, in a coordinate 
basis. Except for trivial cases, this leads to complicated expressions for the Ricci 
tensor that arc undesirable to work with when solving The Einstein field equations. 
We therefore take another approach by noting that we can define another set of basis 
vectors as follows: At a point p £ M, we choose a basis ej and define what is known as 
a left-invariant frame by Lie transporting it around the space in the manifold. That 
is, by the definition of the Lie derivative: 

<» (2-5) 

In general, a frame with basis vectors e,, defines a left- invariant frame, while the 
frame with Killing vectors £j as the basis vectors defines a right-invariant frame. In 
addition, applying the general property: 

£ c . [e i> e fc ]=0 (2.6) 

We see that the basis vectors span a Lie algebra. Therefore, for structure con- 
stants Cy , they obey the general property: 

[e i ,e j ] = C* j e k (2.7) 

We will therefore construct a homogeneous space by first defining a left-invariant 
frame as: 

[e*.ej] = C%B k 

By definition, if oj k is the dual basis to e k , then, necessarily: 

dw* = -ic§w*Ac^ (2.8) 

These invariant one-forms can be used in turn to define a metric (which is also 
left-invariant) as: 

ds 2 = gijuj 1 ® ui j (2.9) 

This is a spatially homogeneous metric on the space. It is interesting to note that 
the Killing vectors of the metric are given by while, the basis vectors ej constitute 
a left-invariant frame. 
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2.2 The Bianchi Classifications 



So far, we have considered the notion of spatial homogeneity in a general sense. We 
obviously would like to apply the aforementioned concepts to general relativity, and 
specifically, cosmology. We begin by assuming that four-dimensional spacetime can be 
foliated with three-dimensional spatial sections. That is, given a spacetime manifold 
M , we have: 

M = Rx S t (2.10) 

Here, M represents the space of time vectors that are moving through the spatial 
slices, and S t is a homogeneous space of dimension three that is unique at each time 
ti. We have attempted to illustrate this concept in Fig. (2.1). 



Figure 1: 3+1 Decomposition into Spatial Slices Orthogonal to the Time Vector 




For a homogeneous space, there are two interesting cases with respect to cosmology. 
The first occurs for dim{Isom(M)} — 6, which means that E t is multiply transitive 
and is maximally transitive. That is, it represents a spatially homogeneous space 
that is translationally and rotationally invariant, or, it is spatially homogeneous and 
isotropic. The models of the present-day universe, the Fricdmann-LcMaitre- Robertson- 
Walker models belong to this classification. With respect to early-universe cosmology, 
we are interested in cases where dim{I som(M)} = 3, which implies that E t is simply 
transitive. That is, we are studying spacetimes that are only spatially homogeneous, 
neglecting the additional assumption of isotropy. The question seems to be then how 
many different classifications of Lie algebras can we have in three dimensions? Such 



8 



classifications are called The Bianchi models, first developed by Luigi Bianchi [5]. One 
typically begins by listing the Bianchi models in terms of their structure constants. In 
general, we can employ the Behr decomposition which enables us to write the structure 
constants as a decomposition of the symmetric and asymmetric parts: 

C k = e l3l n lk + ai (5\5 k 3 - 5]5 k ) (2.11) 

In this equation, aj is a vector belonging to the Lie algebra, for which we have the 
freedom to choose a basis such that a, = adf . 

duj k = --C^w* AoJ j (2.12) 
One then writes the metric tensor as: 

ds 2 = -dt 2 + gij (t)w l ® oj j (2.13) 

This approach is useful when considering the Hamiltonian/Lagrangian formulations 
of general relativity with respect to Quantum gravity, or gravitational chaos. Such 
approaches have been discussed in detail by Ryan and Shepley [31], Hobbill et.al [20], 
Misner [28], Burd [8], and others. Now that we have cemented the geometric structure 
of a spatially homogeneous spacetime, we now apply these concepts to deriving a form 
of The Einstein field equations that will result in a first-order system of ordinary 
differential equations. First, however, we will need a way to describe the presence of 
viscous fluid matter in spacetime. We do this in the next chapter. 

3 The Mathematical Description of a Viscous Fluid 
Matter Source in Spacetime 

So far, we have discussed the topological formalism of a spatially homogeneous space- 
time by introducing the concept of The Bianchi types. This is obviously only one half 
of the story. To develop a general relativistic theory of the early universe, one also 
needs to introduce a matter source that couples to the curvature of spacetime. We are 
particularly concerned with viscous fluids in this thesis. It is believed that the early 
universe contained a matter source that was a high-temperature viscous fluid. [7] [6] It 
is quite plausible that such a fluid had strong properties of turbulent and compressible 
flow. Therefore, in order to account for these effects, it is necessary to include the 
effects of viscosity in any matter source description via the energy-momentum tensor. 

3.1 Energy-Momentum Tensor for a Viscous Fluid 

As per the aforementioned requirements, we need an expression for the energy mo- 
mentum tensor T a ;, for the viscous fluid flow. One is already quite familiar with the 
energy-momentum tensor of a perfect, non-viscous fluid: 

T ab = u a u b {ji+p)-u c u c g ab p (3.1) 
(Where fi represents the total mass-energy density of the fluid) 
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This comes from the idea of considering a pressure-less fluid: 

Too = pu°u° =► T ab = im a u b (3.2) 

Assuming there is a fluid pressure, from elementary fluid mechanics it is well known 
that such a pressure gradient only is defined for the spatial coordinates, so any such 
contribution to the energy-momentum tensor should be proportional to ph ab , where 
h ab is the projection tensor, and is defined as: h ab = g ab — ^j-- We can therefore 
write that: 

T ab = pu a u b + aph ab (a = ±1) (3.3) 
This in turn implies that: 

T ab = u a u b (fx - ap/u c u c ) + apg ab (3.4) 

Letting a — —u c u Cl we find that the appropriate modified form of the energy- 
momentum tensor for a perfect fluid is evidently: 

T ab = (fi + p)u a u b - u c u c g ab p (3.5) 

Letting p, + p = w: 

T ab = wu a u b ~ u c u c g ab p (3.6) 

This is precisely the equation for the energy-momentum tensor for a perfect fluid as 
given in [23]. We proceed with our derivation by first denoting the additional viscous 
contributions by V ab , so that: 

Tab = wu a u b - u c u c g a bP + V ab (3.7) 

To obtain the form of this additional tensor term, we note that from classical fluid 
mechanics, the Euler equation is given as: 

{fiUi) t = -Hik,k (3-8) 

(Where, 11^. ^ is the momentum flux tensor). Also, recall that for a non-viscous 
fluid, one has: 

ILfc = p8 lk + puiu k (3.9) 
We simply add a term that represents the viscous momentum flux, T>' ik : 

n ife = pS lk + pu^k - S- fe = -Sj fe + puiU k (3.10) 

Where 

S ife = -pSik + 24 (3.11) 
is the stress tensor, while, T,' ik is the viscous stress tensor. 
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The general form of the viscous stress tensor can be formed by recalling that vis- 
cosity is formed when the fluid particles move with respect to each other at different 
velocities, so this stress tensor can only depend on spatial components of the fluid 
velocity. We assume that these gradients in the velocity are small, so that the mo- 
mentum tensor only depends on the first derivatives of the velocity in some Taylor 
series expansion. Therefore, T,' ik is some function of the u^fc. In addition, when the 
fluid is in rotation, no internal motions of particles can be occurring, so we consider 
linear combinations of u^k + Wft,i, which clearly vanish for a fluid in rotation with some 
angular velocity, Hi . The most general viscous tensor that can be formed is given by: 

S' ife = r) (u it k + Uk,i - ^S ik ui^ + £5ikU lt i (3-12) 

Here, rj and £ are the coefficients of shear and bulk/second viscosity respectively 
[23] [21]. In the above equation, we note that 5ikUi t i is an expansion rate tensor, and 
(u^k + Uk,i — §<5ifcM/,z) represents the shear rate tensor. Since we would like to gener- 
alize this expression in a general relativistic sense, we replace the partial derivatives 
above with covariant derivatives, and the Kroenecker tensor with a more general metric 
tensor, that is, Sik — > g%k- Therefore, we have that: 

= 77 + u k -i - ^gikui-i^j + igikui-i (3.13) 

As we shall show later when discussing the Raychaudhuri equation, a fluid moving 
through spacetime can be described as a congruence of fluid lines in motion that exhibit 
properties of expansion/collapse, rotation, and shear. Denoting the shear rate tensor 
as a ah-, and the expansion rate tensor as 9 = u a a , the above expression generalizes to: 

V ab = -2 Wab - i9h ab (3.14) 
Substituting these relationships into Eq. (3.7), we finally obtain: 

Tab = (fi + p) U a Ub - u c u c g a bP - 2r)(Xab - £,0h a b (3.15) 

3.2 Equation of Motion for The Fluid 

Now that we have the energy-momentum tensor in hand, we are in a position to 
derive an equation of motion for the fluid. It is quite a remarkable fact that the 
classical equations of motion of fluid mechanics and electromagnetism can be derived 
quite easily from the appropriate energy-momentum tensor. That is, equations of 
motion are typically found from conservation laws. The generalization of conservation 
of energy-momentum in general relativity appears via the divergence-free property of 
the energy-momentum tensor, which itself, is coupled to the divergence-free property 
of the Einstein tensor. 

For notational convenience (which will be clearer in what follows), we will let ir a b — 
2r)<Jab: such that: 

Tab = (P + P)u a U b + gabP - TTafc ~ £,0Kb (3-16) 
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We write this in the mixed form: 

T b = (p+p)u a u b + 5 b a p-7r b a -^h b a 
= ^+p) Ua u b + S b p-n b -^h b a 

= (fi + p)u a u b + S b aP -7T b a -C0(S b a + u a u b ) (3.17) 

Energy conservation requires that this quantity should have a vanishing divergence, 
T a b h = = T b . b . Upon a contraction with u a , we get that: 

u a (( M + p)u a u b + S b p -n b -^e (S b + u a u b )). b = (3.18) 

We will now evaluate this expression term-by-term. For the first term, we get: 



u a [{n + p). b u a u b + (p+p)u a . b u b + (p+p)u a u b b ) 
= u a {{p+p), b u a u b + (fi + p)u a . b u b + (fi + p)u a 9) 
= (-{p + p) fi u b - (n+p)0) 

= (-(i-p-(p+p)8) (3.19) 

For the second term, we have: 



« a (S b , b P + S b a P;t 

'■> 

p (3.20) 



u a (S b aP , b ) 



For the third term, we have: 

= -o ab ^ ab (3.21) 

For the fourth term, we have: 



u a (£6 {S b +u a u b )). b 
= "° (S b + u a u b ) + #. 6 (S b + u a u b ) + £9 (S b + u a u b ) ;b ) 

= -£6 2 (3.22) 

Combining these individual expressions, we obtain: 

(-£ -p-(j* + p)6) +p + <j ab n ab + £6 2 = (3.23) 
Cleaning this up slightly, we have: 
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fi+(ji + p)6- a ab TT ab - £6 2 = (3.24) 
Recall, we had initially defined ir ab — 2r)<r ab , so this equation becomes: 

fi + {^i + p)0 - 2na ab a ab - £,0 2 = (3.25) 
Using the definition a 2 = \a ab a ab => o- ab o ab = 2a 2 , we have finally that: 

fi + (pi + p)6 - Ana 2 - £(9 2 = (3.26) 

We assume that in the early universe, any fluid would behave ultra-rclativistically, 
that is, it would obey an equation of state: 

p = i/i (3.27) 
Substituting this relationship in the previous equation, we get: 

(i + 9 - Ana 2 - £6 2 = (3.28) 

This seemingly simple equation is quite a remarkable result in actuality. It describes 
the evolution of the fluid in our universe. It is important to note, however, that this 
is one equation with three unknown variables. As we shall see later, this evolution 
equation is actually the equation that closes the Einstein field equations constituting 
a complete dynamical system description of the universe and its matter content. 



3.3 A Discussion on Energy Conditions 

Later on, we will attempt to prove the existence of an initial singularity for this viscous 
fluid universe. The arguments presented there will depend strongly on what are known 
as the strong and weak energy conditions. Having just introduced the form of the 
energy momentum tensor, this would be an appropriate opportunity to discuss these 
concepts. Assuming that the universe on a large-scale structure does not have exotic 
matter, that is, matter that does not have a negative energy density, any observer 
moving with four-velocity u a will necessarily measure the energy density (which is a 
scalar quantity) as T ab u a u b . This general condition is typically broken down into two 
specific conditions: The Weak Energy Condition, an The Strong Energy Condition [17] 
[15] . The weak energy condition says that for all u a in a spacetime, 

T ab u a u b > (3.29) 

On the other hand, the strong energy condition says that for all u a in a spacetime, 

(r ab - \Tg ab ^j u a u b > (3.30) 
Using The Einstein field equations, we can write this equation as simply: 

R ab u a u b > (3.31) 
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From a physical point of view, the weak energy condition is a statement about all 
timelike observers measuring a positive energy density. The strong energy condition is 
more restrictive than that in the following sense. To diagonalize the energy-momentum 
tensor, let us choose a frame whose basis consists of the eigenvectors of the energy- 
momentum tensor, and denote the eigenvalues and eigenvectors by /i and pi i = 1, 3. 
Note that it is common to label the p, as the principal pressures. We then see that the 
weak energy condition (WEC) satisfies the condition: 



Furthermore, Eq. (3.31) implies that any spacetime whose corresponding energy- 
momentum tensor satisfies the SEC necessarily has positive curvature. This, in turn, 
says a great deal about the geodesies of this spacetime. In that, for this spacetime 
two neighbouring geodesies initially parallel and separated will eventually converge. 
These concepts will be used later on in the discussion of the necessary criteria for the 
existence of a past singularity. 

4 The Einstein Field Equations 

At the heart of general relativity is the idea that spacetime itself is dynamical. That 
is, in The Einstein field equations is the notion that a matter source induces space- 
time curvature, which in turn, induces motion on the matter source. Since we are 
particularly interested in early-universe cosmology, we will solve these equations for a 
spatially homogeneous, anisotropic, viscous fluid universe. As we alluded to earlier, we 
will specifically work with the Bianchi Type-IV model. The reason for this is largely 
that Bianchi Type-IV models have not been studied in the context of early-universe, 
viscous cosmology. We therefore feel that solutions to this model will provide valu- 
able insights to the fields of early-universe cosmology, and more broadly, cosmological 
dynamical systems. As General Relativity is a coordinate-invariant theory, there are 
often several practical approaches to solving The Einstein field equations. That is, 
there are several appropriate choices for coordinate systems that can simplify the pro- 
cess of obtaining a solution to the field equations. The great detail in describing the 
Bianchi structure above was with great reason. Given the structure of the Bianchi 
models and the general covariance of The Einstein field equations, we will employ the 
powerful method of orthonormal frames. This approached was first used by Ellis and 
MacCallum [14] [42], and in the author's opinion, is one of the most remarkable re- 
sults to come out of general relativistic theory in recent history. In order to use this 
approach, we will first give a covariant description of fluids in spacetime. 



fi>0Afi + pi>0 
Whereas, the strong energy condition (SEC) satisfies the condition: 



(3.32) 




(3.33) 
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4.1 The Orthonormal Frame Approach 

If one considers a fluid described by a family of time-like curves, with 4- velocity, u a , 
the properties of the fluid flow are most appropriately described by the decomposition 
[11]: 

u a ;b = -a a u b + <7 ab + u ab + ^h ab 9 (4.1) 

In this equation, 9 = u a a is a measure of the divergence of the family of time-like 
curves, a a = u a -. b u b is a measure of how much the curves represent non-geodesics, thus, 
can be taken to be the fluid acceleration, and h ab = g ab — u a u b /u c u c is a projection 
tensor. We also have terms that represent the shear and vorticity of the fluid: 



0~ah 



1 , \_l 

2 (^"m;n ^"n;m) ^ 



KK (4-2) 



In addition, we note that since for any velocity field u a u a — ±1, the following 
addition relations hold: 

a c u c = Lo nc u c = a nc u c = h nc u c = (4.4) 

When using orthonormal frames, the metric tensor, g uv takes a very simple form 
for the following reason. Any cosmological model can be described by a coordinate 
system in which we can define a basis {e„}, and its dual basis of differential one- 
forms {w 11 }. As usual, the metric tensor is then defined as g uv — g(e u ,e v ), from 
which the line element is ds 2 — g uv uj u uj v . For a given coordinate chart on the pseudo- 
Riemannian manifold, we can take the basis formed by {e u } as the coordinate basis 
{d/dx 1 }, with the dual basis then playing the role of the coordinate differential one- 
forms, {dx 1 }. If we now assume that the basis constitutes the orthonormal frame of 
this section, then all four basis vectors {e„} are mutually orthonormal, and thus satisfy 
the relationship: g(e u , e v ) = r/ uv = S uv = diag(—l, 1, 1, 1). Hence, g uv is really just the 
Kronecker tensor, and all tensors in this section in the orthonormal frame approach 
have their indices lowered/raised with it. Furthermore, the structure constants are 
now functions, which we denote by a lowercase c to distinguish between the constants 
described previously by an uppercase C. These functions satisfy the commutation 
relation: 

Recall that for an arbitrary basis, for a general (torsion-free) spacetime, we have: 

c a =T a - T a (4 6) 

We can then write the connection coefficients as: 

r a «t) — ^ (Saf> c tm 9uv c av ~ 9v b C a u) (4-7) 
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Building on the discussion of the previous section, we shall assume that the fluid 
is non-tilted. That is, the fluid four-velocity is orthogonal to the spatial hypersurfaces. 
This immediately implies that the fluid is irrotational and geodesic: 

Uuv = u u - v u v = (4.8) 

We can therefore write that: 

@uv l^u-.v "T^Qhuv ~t~ @uv (4-9) 

As we mentioned the four-velocity is orthogonal to the spatial slices, which addi- 
tionally implies: 

0uv = r*„ c\ a = 4 = (4.io) 

As we did with the Behr decomposition, we can write the structure coefficients as: 

c a tb = -9 a b +e a bc W (4.11) 

The vector il c has the special significance of being interpreted as a local angular 
velocity in the rest frame of an observer with spatial frame {e a }. It is conventionally 
defined as: 

W = \e ahcd u b e c ■ e d (4.12) 

It can be shown that that the structure coefficients are spatial, and therefore, must 
correspond to The Bianchi models. We can therefore write that: 

4 = e ijin lk + ai {5\5 k - 5$) (4.13) 

Since each structure coefficient is constant along each orbit of transitivity, n ab and 
dj are functions of time only. Even though the spatial frame described by the basis 
{e u } is a left-invariant set on the spatial slices, the difference in the orthonormal frame 
approach is that the structure coefficients are time-dependent functions. Before, the 
structure coefficients were always taken to be constant in the true sense of the word. 
We can find evolution equations for n ab and at by first realizing that the Jacobi identity 
holds for all vectors. The Jacobi identity as applied to the set of vectors (U,e a7 e b ) 
gives: 



[U, [e a , e 6 ]] + [e a , [e 6 , U]] + [e b , [U, e a ]] = 

=> (U« 6 ) + « b + « t + cl u c u ta )e v = (4.14) 

In this derivation, U is a gauge vector that will typically be set to dt- Note that since 
we are assuming that the fluid is non-tilted, necessarily 6 UV — => c\ a — c^ fc = 0. 
We therefore obtain the identity: 

dt(c k ab ) + c k td c d ab + c k ad c d b t + c k bd4a = (4.15) 
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One can show that upon applying Jacobi's identity to the three spatial vectors, we 
get the eigenvalue equation [15]: 

n ij a,i = (4.16) 

Combining the previous equations and taking the trace of Eq. (4.15) gives the 
evolution equation for a,i as: 



The trace-free part of Eq. (4.15) gives the evolution equation for n a b as: 



(4.17) 



n a b 



-0n„ 







(4.18) 



The important thing to keep in mind is that we would like the structure constants 
as defined in Eq. (4.13) to correspond to a Lie algebra. In order for this to be the 
case, Eq. (4.16) must hold. Bianchi models are typically classified into two categories: 
Class A models and Class B models. Class A models are all Bianchi models such that 
di = for which the latter is satisfied. For class B models, a% must be an eigenvector 
of n lJ with zero eigenvalue. It is important to note that we always take n y to be a 
symmetric matrix and as such we can diagonalize it using a orientation of our choice 
for the spatial frame. It has been conventional to assume: 



riij = diag(n±,ri2,ns), , a 1 = (0,0, a) (4-19) 

The Jacobi identity immediately implies that n 33 a = n 3 a = 0. It is a fact of linear 
algebra that eigenvalues of a matrix are invariant under a conjugation operation with 
respect to rotations. One can then classify the different Bianchi models by the signs of 
the eigenvalues «n, «22i "-33 = ni,n2,ris and a. Since, we are dealing with the Bianchi 
Type IV, we note that: 

a 1 = a5l > 0, and m > 0, n 2 = n 3 = 0. (4.20) 

(The signs of the eigenvalues rii and a for the other Bianchi types can be found in 
[42] [15]. 



4.2 The Dynamical Equations as Implied by The Einstein Field 
Equations 

In this chapter, we have so far explained the orthonormal frame approach, derived 
the evolution equations for n,j and a % , and described the classification of the Bianchi 
Type-IV model. We will now derive the dynamical equations which are implied by 
The Einstein Field equations in the orthonormal frame approach. To give a sense of 
how remarkable this approach is to solving the field equations, one should first look at 
them written in terms of the metric tensor alone [39] . 
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2~ab 



d 2 ~g 



dx u dx 



V ~av~bu , ± l7,a.b~ ~Xv ~pu\ , ~ ~vp~al~bu 

~ s,„ + ^ (5 giu9, p g% ) + gi u g g, v g, P 



~bv~up ~bl~ ~av~up 



+ 



g 9uvg, p gf - g gu V g, p g^ 

~ab~lu\ /o 7. 7, 7, 7. \ TjVt-ps 



(2g av ~g bu - g ab g lu ) (2& P &t - ~g ps g v t) gfg% 



(4.21) 



This system of partial differential equations represents the vacuum field equations 
in the deDonder gauge (g a b = 0), where g uv = \f-~gg uv ■ This is a very complicated set 
of ten coupled highly nonlinear partial differential equations. As we shall show below, 
the field equations lead to a system of first-order ordinary differential equations in the 
orthonormal frame approach. 

4.3 The Raychaudhuri Equation 

We begin by recalling that one can consider a fluid evolution as described by a family 
of time- like curves, with 4- velocity, u a . We had from before: 

u a ;b = -a a u b + <J ab + u ab + ^h ab 9 

Analyzing the equation of the divergence of the fluid flow, we are interested in how 
this divergence evolves in time: 



d9 



= (u? a ). b u b = u^ b u b 



(4.22) 



Recalling the definition of the Riemann curvature tensor as determine whether the 
second covariant derivatives of some vector field in the manifold do not commute: 



U a ;d:b ~ 


Ua;b;d 


— R a db u c 


=> <d;6 


u ;b;d 


= Rdb u c 


=> u °d;b 




= -R%uc 


=> u -d-b 


u ;b;d 


= ~R-cdb u 



Performing a contraction on the indices a and d, we finally get: 



;a;b 



■b-a 



R bc U C 



Substituting this equation into Eq. (4.23), we obtain: 

^ = u%. a u b - R bc u c u b 
Using the fact that: (^u a b u b ^j = u a b . a u b + u? b u b a , we have: 



(4.23) 



(4.24) 



(4.25) 
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d9 



= a a a - R hc u c u - u a . h u 



Since h a bh ab = 3,er 2 = \(J ab a a b, lu 2 = ^Lu ab uj ab , we now have: 

u a:b u b -- a = 2o 2 - 2u 2 + U 2 + hh ab a ab 
In addition, we also have that: 



h ab u a , b = -a a u b h ab + h ab a ab + h ab Lo ab + -6h ab h ab 



Since: h ab a a b — 0, and that: 



u a , b u b ' a = 2a 2 - 2uj 2 + -0 2 



We now finally get the Raychaudhuri equation for 9: 



d£_ 
~dt 



-R bc u b u c + o° + 2lu 2 - 2a 2 - -9 2 

-R bc u b u c + (o° - T a an a n ) + 2uj 2 - 2<i 2 - ^9 2 



We can make a further simplification using the Einstein field equations 
for the Ricci tensor term in the above equation: 



R a b — -^Rg a b 

^ R-2R 
=> R 

Rab + 2 KT 9ab 



nT ab 

kT 
-kT 



Rab — 



K[T a b- -^Tg a b 



We compute the trace of the energy-momentum tensor as follows: 

T = T a a = g ab T ab = g ab (wu a u b - u c u c g ab p - 2 Wab - ^h ab ) 

= g ab (wu a ub) - g ab (u c u c 9 abP) - 2 V g ab a ab - ig ab (9h ab ) 

= u a u a (w - P S a a ) - 2 V g ab a ab - Zg ab (9h ab ) 

= u a u a (w - 4p) - 2 V g ab a ab - Hg ab (9h ab ) 

= u a u a {n - 3 P ) - 2 V g ab o- ab - Hg ab (9h ab ) 

We can simplify things slightly by writing: 



0~ab 



1 1 



Kt'K = u {a:b) + a (au b ) ~ -^0h ab 
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We then have: 



g ah <j a b = g ah \ U( a -,b) + a ( a u b ) - ^QKb 

= g ah u { a:, b )+9 ab a Uub) - l -e 9 ah h ab 

= U-a + 0-6 

= (4.34) 



Here we have used the fact that: 



ab h ab ( „ U a u b 

g hab = g g a b 



u"-u c 



a U C U C 



We then have: 



(4.35) 



g ab (6h ab ) = 6g ab h ab = 39 (4.36) 
Therefore, the trace of the energy-momentum tensor T takes the form: 



(4.37) 



T = u a u a ( f i-3p)-2r,g ab a ab -^g ab {0h ab ) 
= u a u a (n - 3p) - 3£9 

From our original definition of T ab , it follows that: 

T ab u a u b = (wu a Ub - u c u c g a bP - 2Wab - £0h ab ) u a u h 

= wu a ubu a u b — u c u c g a bu a u b p — 2r]a a bU a u b — £,6h ab u a u b 

= w {u a u a f -p- 2 Wab u a u b - £,6h ab u a u b 

= w - p- 2ri<j ah u a u b - £9h ab u a u b 

— w — p 

= H+P-P 

= fi (4.38) 
Combining the results we obtained above, we obtain: 
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RabU a U b = 



nu a u° ( T ab - -Tg ab 



nu a u b T ab - - K Tg ab u a u b 

- 7} K [u a u a {p, - 3p) - 3£,9] u a u a 

- X -n [(ji - 3p) - 3^u a u a ] 

13 3 
up - -Kfi + -up + -n^6u a u a 

1 3 

-k (p + 3p) + -n£6u a u a 



Since we committed to the signature (—1, +1, +1, +1), u a u a = — 1, and: 



(4.39) 



R ab u a u b = -«(/x + 3p)- -k£0 
Therefore, we obtain 27ie Raychaudhuri equation for our viscous fluid: 



(4.40) 



§■ = -R bc u b u c + ° + 2u? - 2a 2 - \e 2 

-\k {p + 3p) + §/s£0 + a? a + 2to 2 - 2a 2 - \o 2 



Since we are assuming the fluid is geodesic and irrotational, we obtain: 
d9 1 , „ N 3 _ „ o 1 



dt 



- K ^ + 3p) + ^6-2a 2 --e 2 



Since p = ip, we finally obtain: 



p 2 3 



(4.41) 



(4.42) 



(4.43) 



Note that we have taken some care to derive The Raychaudhuri equation by going 
through each step carefully since such a derivation for a bulk viscous fluid is not readily 
available in the literature. 

4.4 The Generalized Priedmann Equation 

Exploiting the assumption of a non-tilted cosmology again, we state the additional fact 
that in this case: 



"a;6 — ab — K ab 



(4.44) 



21 



Here, we have denoted the extrinsic curvature of the spatial slices by the extrinsic 
curvature tensor, K ab . Recall from standard differential geometry, we have that for 
any three-dimensional spatial slice: 

WR = ^R + K 2 -K a ^K aP -2^R a0 u a u^ 
=>■ kT = - (3) R- K 2 + K aP K a p 



2^R a0 u a u fi 



/,T = -WR- K 2 + K al3 K a/3 + 2KT a pu a u fl + kT 







= - (3) i?- K 2 +K ali K aP + 2KT af} u a u' 3 



= \^R-K^K aP + K 2 ) 



(4.45) 



If we now take once again the decomposition equation and set a a = cj a j = 0, we 
obtain: 



Ua;/3 — 0a/3 = K a p = (J a p + -h a /30 



In addition we note that: 



K 



1 



K afi K a ? 



(4.46) 



(4.47) 



Substituting Eq. (4.38) into Eq. (4.45), we obtain: 



Kfl 



= l -^R-K^K af) 



K- 



(4.48) 



Using the dchnition in Eq. (4.47), one obtains the generalized Friedmann equation: 

(4.49) 



II 1 (3) 

1 Q2 = 2 aaha 2 R + K ^ 



4.5 The Shear Propagation Equations 

The last set of dynamical equations implied by The Einstein Field equations are the 
shear propagation equations. They essentially describe the evolution of the anisotropy 
in a cosmological model as a function of time. The lengthy derivation of these equations 
has been done in much of the literature on cosmological dynamical systems [42] [15] [29]. 
The idea behind the derivation is similar to that of the derivation of The Raychaudhuri 
equation. The shear propagation equations are: 



cr'ab + 0o- ab - <r d a e bcd n c - ate acd n c + (3) R a b - ~4& )jR = 2 V^ ab 



(4.50) 
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4.6 The Ricci Curvature and Constraint Equations 



The greatest achievement of the orthonormal frame approach is that one no longer 
needs to express the components of the Ricci tensor in terms of the coordinate basis 
functions. The Ricci tensor is expressed entirely in terms of a % and n a b. As discussed 
by Ellis and MacCallum [14] and Gr0n and Hervik [15], the Ricci tensor takes the form: 

{3} Rab = -£ c a d n bc a d - e c b d n ac a d + 2n ad n d - nn ab - S ab ^2a 2 + n cd n cd - ^n 2 ^j (4.51) 

Contracting this equation, we obtain the Ricci scalar: 

^R ^ R a = - [&a 2 + n cd n cd - ^ (4.52) 

One can also show that the off-diagonal components of the four-dimensional Ricci 
tensor yield a constraint equation: 

3aV fca - e abc n cd a b d = (4.53) 
It will however be more convenient to write this equation out in component form: 

3acr 33 + (n 11 - n 22 ) a 2 i = 
3a<T 3 i + n 22 a 32 = 

3acr 32 - n n cr 31 = (4.54) 

We have now discussed and derived all the necessary formalism of a spatially ho- 
mogeneous spacetime that is additionally anisotropic. We will now finally apply these 
results to The Bianchi Type IV metric in the next chapter. 



5 The Bianchi Type IV Dynamical Equations 

In this chapter, we will derive the dynamical equations associated with the Bianchi 
Type IV algebra, and study their solution. For the Bianchi Type IV geometry, we take 
nn = N, and 7122 = w 33 = 0, where N > 0. The constraint relations (Eq. 4.54) then 
lead to: 



3a<7 33 + NC21 = (5.1) 

3acr 3 i = (5.2) 

3ao- 32 - 031 = (5.3) 

We can see that tr 3 2 = 0. 

Since the shear tensor is taken to be symmetric and traceless, in three dimensions, 

it only has two independent components. It is also clear that the only non-zero off- 
diagonal component is 021 = &12- 
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Therefore, we have that, for a = b: 

a ah = (a + + V3cr_ , cr + - V3cr_ , -2<r + ^ (5.4) 

We have also one non-zero, independent, off-diagonal, component for the shear 
tensor: 

-3acr 33 6aer + 

*21 = <?12 = — ~— = — (5.5) 

Using Eqs. (4.17) and (4.18) we obtain: 

n 1 = n 2 = o (5.6) 

a + ^Oa + cr 33 a = 

a + ]-0a - 2<7 + a = (5.7) 



N- ^9N -2N (a + + V3a-^j = (5.8) 

fl 3 = -6aa + (5.9) 

We also compute the components of the three-dimensional Ricci tensor as follows: 

{3) Rab = -e c a d n bc a d - e c b d n ac a d + 2n ad n b - nn ab - h ab ^2a 2 + n cd n cd - ^n' 
= -fl 3 "(ii(fl) ~ 4 3fl ai(a) + 2n al n\ - nn ab - S ab [2a 2 + nun 11 - -n 2 



-e a 3 n bl (a) - e\ 3 n al (a) + 2n al n\ - Nn ab - S ab ^2a 2 + ^N 2 ^j (5.10) 



Now, if a = 1,6 = 1: 



(3) i?n = -N 2 - 2a 2 
2 



(5-11) 



If a = 1,6 = 2: 



If a = 1,6 = 3: 
If a = 2, 6 = 2: 



( VR 12 = (3) R 2 i = -efn n a 

= Na (5.12) 

^R 13 =W R u = (5.13) 
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(3) i?22 = - ^2a 2 + V 

If a = 2, b = 3: 

(3) i?23 - (3) i?32 = 

If a = 3, 6 = 3: 

(3) #33 = - (2a 2 + X -N 2 
The Ricci scalar is computed from Eq. (151) to be: 

R =W Rl = -6a 2 - -N 2 
Using Eq. (4.50), we derive the shear propagation equations as follows: 

<Jab + 8o ab -<j d a e bcd n c -a b d e acd n c +^ R ab - ls<$ R = 2 V >ia ab 
v'ab + 0a ab - a d a e b3d n 3 - ate a3d Cl 3 +^ R ab - ^ b } R = 2 V K<j ab 
^ <y'ab + 0<J ab - ^ 3 (<J^ b3d + 4e a3d ) +^ R ab - -5 a 3 b ] R = 2 V K<J ab 

For the diagonal equations (a = b) wc obtain: 

/ \ 27V 2 / \ 

cr'+ + V3a- + 6[a+ + Via-) - 12a 2 a\ + = 2r)K (a + + V3er_J 

<T+ - V3a- +e(a + - V3cr_) + 72a 2 cr 2 - ^- = 2t)n (a + - V3<7_) 



N 2 

<t+ + a +6 + —— = 2r]na + 


As for the off-diagonal equations, we obtain: 



6aa + + 6aa + — h 69aa + + 12v / 3acr+o'_ + N 2 a = 12-qnaa + 

Recall, we had that: 

a = — a ( — 9 — 2a + 
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N = ^6N + 2N (o+ + V3cr_) 
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Substituting these results into Eq. (5.22), we have upon simplification: 

iV 2 + 6er+ + 2a+6 = 12k W+ (5.23) 
If we now subtract Eq. (5.19) from Eq. (5.20), we obtain: 

N 4 - lAAa 2 a 2 + + 2V3N 2 [ct_ + a- (9 - 2ktj)} = (5.24) 
We can also add Eq. (5.21) and Eq. (5.23) to obtain: 

-5iV 2 

5cr+ - a+6 = -10nr]<7 + (5.25) 

6 

Equations (5.24) and (5.25) therefore constitute a coupled system of first-order, 
ordinary differential equations for the shear variables a + , cr_ . 

We will additionally write the generalized Friedmann equation in the following 
convenient form: 

= i(9 2 - 3a 2 - jN 2 - K[i (5.26) 

Using this relationship, we will write the dynamical equation for the fluid in the 
following form: 

£=-12o 2 rj-4s W -»?JV 2 -^ + ^+0 2 £ (5.27) 
We will also write The Raychaudhuri equation in the following form: 

9 = - l -K( j i + ip)+ i -^9-2a 2 -U 2 

= i (I2a 2 + 2/t/i + N 2 - 26 2 + 3/c0£) (5.28) 

Therefore, our full cosmological system for the Bianchi IV viscous fluid universe is 
governed by a system of six, coupled, ordinary differential equations in six unknowns, 
with the two viscosity coefficients, i] and £: 
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2k^i + N 2 



V = 



-\2c?r\ - iK^rj - r)N 2 



N 2 



24V3a 2 a 2 _ 



cr+ = - 



2^3 N 2 
N 2 (7+6 
~6 5~~ 



28 2 + 3k9£) 

~ ~3~ + "IT 

a -6 + 2n<j-r\ 



2nr]a+ 



-a I - - 2(j^ 



N 



N 



^(7+ + V3(T- 



(5.29) 
(5.30) 
(5.31) 

(5.32) 
(5.33) 

(5.34) 



Clearly, this system of equations has no exact solution, so numerical methods must 
be applied. However, even for a six-dimensional system, such as this one, numerical 
algorithms can be difficult to employ to obtain any relevant solutions. Therefore, we 
will write these equations in dimensionless form, by writing them using the expansion- 
normalized variables [42]. 

5.1 Dynamical Equations in Expansion-Normalized Variables 

The basic idea of the expansion-normalized variables is that the set of differential 
equations in Eqs. (5.29-5.34) is of the form: 



dt 



= /(x) 



(5.35) 



where x = (8, /i, cr_ , cr+ , a, N) £ R 6 . We will define 6 = 3H, where H is denoted 
the Hubble scalar and is defined as: 



H = - 

s 



(5.36) 



where s is a cosmological length-scale function. It is also convenient to define the 
cosmological deceleration parameter q: 



ss 

1 = -TT 



Clearly, we have the relationship: 

H = -{l + q)H 2 

It is also necessary to introduce a dimensionless time variable r as: 

s = sqc/ 



(5.37) 

(5.38) 
(5.39) 
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From these equations, one can show that: 

£-i*£-*-~(i + ,)H <"») 

Substituting = 3i7 into the Raychaudhuri equation, we obtain: 

3H = ~K{ji + 3p) + ^«£ff - 2cr 2 - 3i? 2 
^>H = ~K(p + 3p) + «£ff - ^o- 2 - ff 2 

^>H = -\kh + k£H- \g 2 -H 2 (5.41) 

In addition, we define: a density parameter, 

0=^ (5.42) 
a shear parameter, which compares the rate of shear with the rate of expansion, 

* - j£ (5-43) 

and a curvature parameter: 



Oil" 

Note that, if the dynamical system has an evolution equation for fl, then this 
parameter is not explicitly needed in the system of differential equations. We have, 
however, just included it for completeness. 

In addition, comparing Eqs. (5.41) and (5.38) we have the relationship: 



(l + q)H 2 = -\kh + k£H- \o 2 -H 2 
(l + q)H 2 = - + |<7 2 + H 2 

M 1 < 2(j2 , 

=►(1 + 9) = sm^-H-3HS + 1 
1 < 2cr 2 



^q - 3i/2 « M H ^ 

^>q = kQ - - 2S 2 (5.45) 
H 

We additionally assume the bulk viscosity obeys the equation of state: 

| = 3£o (5.46) 
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Here, £o is a positive parameter. This equation of state essentially says that as the 
universe expands the bulk viscosity coefficient asymptotically vanishes describing the 
present-day, perfect fluid universe. 
Therefore, we get: 



q = kQ - 3/c£o - 2S 2 (5.47) 

Why this approach is so powerful in terms of analyzing the dynamics of our system 
is for the following reasons. Eq. (5.30) which is the evolution equation for the energy 
density, is now replaced (in the dynamical system) by 51. Based on physical arguments 
of physical matter having a nonnegative energy density, this quantity is taken to be 
positive. Therefore, the solution space of the system of differential equations has the 
constraint that it is effectively bounded by f2 > 0. Keeping these concepts in mind, 
we will now write Eqs. (5.31) - (5.34) in dimensionless form. (Below, we define the 
additional dimensionless quantities: n = N/H, A = a/H, S± = <j±/H): 

Beginning with Eq. (5.31), we get: 



n 2 H 2 r A 2 U 2 Y?,H 2 

E_H = — + 24V3 —t 3£_ff 2 + 2 K S_Fr? 

2V3 n 2 H 2 

2tt2 A 2 B 2 Y? H 2 

+ ffE_ = ^+24\/3 ^4 3E_iJ 2 + 2KS_iJ7? 

2^3 n 2 H 2 

2 tt2 A 2 B 2 Y? H 2 

=^£_i? = ^+24V^ ^-4 3E_F 2 + 2kZ-Hti - 

2^3 n 2 H 2 

n 2 H n A 2 H 2 E 2 H 2 H 

=>XL = = + 24^3 ± 3S_iJ + 2kS_77- — E_ 

2\/3 n 2 H 3 H 

n 2 H r A 2 Y?\H H 

2%/3 n 2 H 

(5.48) 

Using the chain rule, we obtain: 

d£_ dE_ dt ■ 1 , 

Therefore, we have that: 



= 24^3^4 -3£_-2k£_-^--^E_ 

2y/3 n 2 H H 2 



2 A2^2 



A 2 Z 2 n H 

+ QV 1,,V ' 



2v2 



-A 2 I] 



24^^-4^ - 3E_ - 2kE_-^ + (1 + g)E 



2\/3 n 2 H 

A 2 Y? 

_ , 24\/3 =± _ 3S _ _ 6kS_t; + (1 + g)S_ (5.50) 

2V3 " 2 



As is done by Saha and Rikhvitsky, [32], we have set in the above equation, r//H = 3rjo, 
where t)q is a nonnegative parameter. 
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Continuing with Eq. (5.32), we have: 



E+i? = 



E+7J 



=^E' + = 
=►2' = 



- l -(n 2 H 2 )-^ + H 2 ) + 2n^ + H 

-l(n 2 H 2 ) - l(^+H 2 ) + 2k V Z+H - £TE+ 
6 5 



1 



H 



-~{n 2 H) - -E+if + 2 K nE+ - -E+ 

-in 2 - jjs + + 2k-|e + + (1 + 9 )E+ 

-in 2 - ^E+ + 6Kno£+ + (1 + g)E + 
5 



Going on to Eq. (5.33), we have: 



(5.51) 



AH = -AH (H -2Y.+H) 
AH + HA = -AH 2 + 2E+ Aff 2 



H 



AH + 2E+Aff A 

+ H 



^>A = 

=> A' = -A + 2Y.+A+ (1 + q)A 
Finally going to Eq. (5.34), we have: 

nH = nH H + 2 (z+H + V3E_.ff) 
^>hH = nH H + 2 (e + H + V3Y,-H 

H + 2 (p+H + V3S_ff) 

H + 2 (z+H + V3E_ff 
2n (e+ + V3E_) +{l + q)n 



n = 



n 

n 
H 



-Hn 

H 

H 71 
H 



(5.52) 



(5.53) 



Therefore, the non-dimensional dynamical system of equations takes the form: 
- 2 -A 2 Y?, 



E' = 



K = 



2^3 



2AV3 ^- - 3E_ - 6/cE_7/o + (1 + ?)E. 



1 3 

-n 2 - -E + + 6k?7oE + + (1 + g)E+ 



A' = -.4 + 2E+.4 + (1 + <?),4 

n' = n + 2n(E + + V3E_) + (1 + gr)r 



(5.54) 

(5.55) 
(5.56) 
(5.57) 



One can see the advantage of this approach immediately. The solution space of the 
system has been reduced from R 6 to R 5 . (The remaining evolution equation for the 
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energy density is derived below). Raychaudhuri's equation has been decoupled from 
the system and is not necessary to solve to determine the dynamics of the cosmological 
system. It is important to note that (E + , E_) describe the anisotropy in the Hubble 
flow, and (A, n) describe the spatial curvature of each spatial slice of the spatially 
homogeneous spacetime. 

We now use Eq. (5.43) to calculate S 2 as follows: 



1 



s2 = (< 7 abv ab ) 



6H 2 
1 

6H 2 



1 (a 2 + 3a 2 + a\ + 3a 2 + 4<7 2 ) 



6H 2 
1 



= + (5.58) 
In addition, K is computed from Eq.(5.44) as follows: 



K 



' m 2 

1 Ua 2 + In 2 



6H 2 V 2 



= A 2 + ^ (5.59) 

The O term in the above equations is governed by the non-dimensional version of 
Eq. (5.30), which we obtain as follows: 



2 



i}(3H 2 ) + GflHH = -12H 2 A 2 r t -l2Kr 1 H 2 n-rin 2 H 2 -12H 3 n + 12riH 2 + 9H 2 £ 

^Cl(3H 2 ) = -12H 2 A 2 7]-l2Kf 1 H 2 n-T 1 n 2 H 2 -12H 3 n + 12r]H 2 + 9H 2 ^-6flHH 

n 2 

= -AA 2 r] - 4/tr/O -r)—- 4Hn + 4?? + 3£ + 20(1 + q)H 



=>0' = -12A 2 tj - 12kO??o - mn 2 - 40 + 12% + 9£ + 20(1 + q) (5.60) 
Therefore, the cosmological dynamical system takes the form: 
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£'_ = T + j± - 3S_ - 6 K £_r/ + (1 + g)£_ 

2 v 3 ^ 

£' + = -in 2 -^E + +6« % S + + (l + g )S + 

A' = -A + 2Y l+ A+{l + q)A 

n' = n + 2n(s + + \/3£-) +(l + q)n 

n' = r/ (-12A 2 - 12/tO-n 2 ) -40 + 12% + 9£o + 20(1 + q) 

(5.61) 

where q = k0 - 3< - 2£ 2 , S 2 = E 2 . + E 2 , # = ^ 2 + and > 0. 

6 Solutions To The Bianchi Type IV Dynamical Equa- 
tions 

In this chapter, we present several numerical solutions to the system of equations 
(5.61). Clearly, this system, which is nonlinear and coupled has no exact solutions. Our 
goal however is to discover values for £o and 770 such that the system asymptotically 
isotropizes, that is, S± — > as t — > 0. Such physical characteristics resemble an 
anistropic early-universe that over time isotropizcd to the present-day universe. Note 
that in the diagrams below, asterisks indicate the initial conditions for the numerical 
experiments. 
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6.1 Plots of Numerical Experiments 

Figure 2: Phase plot of the anisotropy in the Hubble flow for £o = 0.1, rjo = 0.1. 
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Figure 3: Three-dimensional phase plot of the anisotropy in the Hubble flow and spatial 
curvature for £o = 0.1, rjo = 0.1. 
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Figure 4: Phase plot of the anisotropy and spatial curvature for £o = 0.1, rjo = 0.1. 




Figure 5: Plots of the Energy Density and Anisotropy variables as functions of time, for 
£o = 0.1, t/o = 0.1. One can see that Q — > 0, £_ — > 0, but S + -/> 0, so the model only 
partially isotropizes. 




Figure 6: Phase plot of the anisotropy in the Hubble flow for £o = 0.1, r/o = 1. 
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Figure 7: Three-dimensional phase plot of the anisotropy in the Hubble flow and spatial 
curvature for £o = 0.1, 770 = 1. 
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Figure 9: Plots of the Energy Density and Anisotropy variables as functions of time, for 
£o = 0.1, r]o = 1. One can see that — >• 0, £_ — > 0, but S + -f+ 0, so the model only 
partially isotropizes. 
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Figure 10: Phase plot of the anisotropy in the Hubble flow for £o = 0.1, r/o = 10. 




41 



Figure 11: Three-dimensional phase plot of the anisotropy in the Hubble flow and spatial 
curvature for £o = 0.1, 770 = 10. 
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Figure 12: Phase plot of the anisotropy and spatial curvature for £o = 0.1, ??o = 10. 
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Figure 13: Plots of the Energy Density and Anisotropy variables as functions of time, for 
£o = 0.1, r]Q = 10. One can see that Q — > 0, but S_ 0, and S + -/$■ 0, so the model does 
not isotropize at all in this case. 
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Figure 14: Phase plot of the anisotropy in the Hubble flow for £o = 0.1, r]o = 100. 
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Figure 15: Three-dimensional phase plot of the anisotropy in the Hubble flow and spatial 
curvature for £o = 0.1, rjo = 100. 
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Figure 16: Phase plot of the anisotropy and spatial curvature for £o = 0.1, rjo = 100. 




Figure 17: Plots of the Energy Density and Anisotropy variables as functions of time, for 
£o = 0.1, r/o = 100. One can see that O — > 0, but £_ ft 0, and £+ ft 0, so the model does 
not isotropize at all in this case. 
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Figure 18: Phase plot of the anisotropy in the Hubble flow for £o = 1, t/o = 0-1- 
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Figure 19: Three-dimensional phase plot of the anisotropy in the Hubble flow and spatial 
curvature for £o = 1 ; Vo = 0.1. 
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Figure 21: Plots of the Energy Density and Anisotropy variables as functions of time, for 
£o = 1, 770 = 0.1. One can see that O — > 0, S_ — > and £+ — > 0, so the model does in fact 
isotropize. 
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Figure 22: Phase plot of the anisotropy in the Hubble flow for £o = 10, r]o = 0.1. 




Figure 23: Three-dimensional phase plot of the anisotropy in the Hubble flow and spatial 
curvature for £o = 10, 770 = 0.1. 
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Figure 24: Phase plot of the anisotropy and spatial curvature for £o = 10, 770 = 0.1. 




Figure 25: Plots of the Energy Density and Anisotropy variables as functions of time, for 
£o = 10, 7/0 = 0.1. One can see that Cl — > 0, S_ — > and E + — > 0, so the model does in 
fact isotropize. 




Figure 26: Phase plot of the anisotropy in the Hubble flow for £o = 1, i]o = 0. 
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Figure 27: Three-dimensional phase plot of the anisotropy in the Hubble flow and spatial 
curvature for £o = 1, 77Q = 0. 




5* 



Figure 28: Phase plot of the anisotropy and spatial curvature for £o = 1> Vo = 0- 




Figure 29: Plots of the Energy Density and Anisotropy variables as functions of time, for 
£o = 1, T]o = 0. One can see that O — > 0, £_ — > and E + — > 0, so the model does in fact 
isotropize. 
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Figure 30: Phase plot of the anisotropy in the Hubble flow for £o = 100, 770 = 0. 
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Figure 31: Three-dimensional phase plot of the anisotropy in the Hubble flow and spatial 
curvature for £o = 100, r/o = 0. 
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Figure 32: Phase plot of the anisotropy and spatial curvature for £o = 100, 770 = 0. 
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Figure 33: Plots of the Energy Density and Anisotropy variables as functions of time, for 
£o = 100, 770 = 0. One can see that Q — > 0, £_ — > and S + — > 0, so the model does in fact 
isotropize. 
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6.2 Discussion of Results 

First looking at Figs. (6.1-6.16), one can see that the cosmological model does not 
isotropize in all directions. In addition, from the phase plots of the dynamic variables, 
one can see that the dynamics of the model resemble that of plane-wave solutions. 
Such behaviour was noted by Hewitt and Wainwright [42] in their investigation of 
the dynamics of a non-tilted perfect-fluid Bianchi Type-IV model. In addition, the 
asymptotic limit of anisotropy of a Bianchi Type IV perfect fluid model was noted by 
Hervik, van den Hoogen and Coley [19]. It is interesting to note that as the shear 
viscosity coefficient 770 increased, the anisotropy in the cosmological model increased 
as well. However, looking at Figs. (6.17 - 6.24) where we introduced a bulk viscosity 
coefficient £0 that was significantly larger than the shear viscosity coefficient caused 
the model to isotropize. What is additionally interesting is that as the bulk viscosity 
was increased, the relationship between the dynamical variables became increasingly 
linear, which hints at the notion that the dynamics of the transitionary period from 
the early-universe cosmology to the present-day cosmological models become simpler. 
Perhaps the most interesting point of all is that while the shear viscosity coefficient 
was set to zero in Figs. (6.25 - 6.32), with the bulk viscosity coefficient becoming larger 
and larger, the model still isotropized. This provides evidence that the presence of a 
significant bulk viscous pressure rather than a significant shear viscosity is necessary 
for such a cosmological model to isotropize. Similar results for other Bianchi models 
have been reported by van den Hoogen and Coley [40], Singh and Kale [35], and 
Pradhan, Rai and Singh [30]. It is also important to note that no cosmological model 
that is spatially homogeneous and anisotropic and necessarily has a three-dimensional 
isometry group can evolve to a model with a six-dimensional isometry group. When 
one speaks of a model isotropizing, £± — ¥ only asymptotically. However, any Bianchi 
model can become arbitrary close to that of the FLRW models as a consequence of 
asymptotic isotropization. It is also a strong implication that an isotropic universe 
necessarily implies the existence of an isotropic fluid. Therefore, one must conclude 
that the nature of our present-day spatially homogeneous and isotropic universe as 
described by the FLRW models is actually quite special indeed! 

7 On The Penrose-Hawking Singularity Theorem 

Up to this point, we have demonstrated that a Bianchi Type-IV model isotropizes with 
the significant presence of bulk viscous pressure. One can interpret our work so far 
as evolving an early-universe cosmological model towards the future. Now, we wish to 
analyze the reverse question by considering the physical past of a viscous model. In 
particular, we are interested to know whether such a universe could have emerged from 
a physical singularity point. 

We begin by considering The Raychaudhuri equation as given in Eq. (4.42): 

= ~k (/i + 3p) + |/s# - 2a 2 - V 
If we assume the strong energy condition (Eq. (3.31)): 
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R ab u a u b > 
=» \n (/i + 3p) - > (Eq. (4.39)) 

-l K ^ + 3p) + ^e < (7.1) 



Here, we assume the equation of state: 



Hence, Eq. (7.1) becomes: 



j,CeR+ (7.2) 



^K(/i + 3p) + |/tC<0 (7.3) 



Since <r 2 > 0, we have from The Raychaudhuri equation: 

+ ^/s (/i + 3p) - \kC+^9 2 = -2a 2 <0 (7.4) 
Using Eq. (7.3), this becomes: 

e + ± K (n + 3p)-^KC + ^e 2 < o 

^>0 + ^6 2 < -^k(h + 3 P ) + ^kC < o 



+ -e 2 < o 



< -lo 2 (7.5) 



We first proceed assuming a strict inequality in Eq. (7.5). That is: 



1 

d fl\ 1 
> - 



|2 



dt \9 J ~ 3 

*m - H (76) 

We have assumed the initial condition: 9i = 0(0). In addition, we also assume that 
t < 0. If we have an expanding universe, then, necessarily, the congruence of geodesies 
is expanding at t = 0, which implies that 0, > 0. This in turn implies that: 

VteM,3 t sing such that: 1 = (7.7) 

vitsing ) 
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That is, the function l/0(t) must have passed through zero at some time, t a i ng . It 
then follows that: 

6{t sing ) = oo (7.8) 



9(tsin g ^) 

This tells us that there was a physical singularity at the time t s i ng , where the 
expansion scalar was infinite. 

In the case of a strict equality in Eq. (7.6), we can be slightly more precise by 
employing the intermediate value theorem. We begin by defining the function: 

We wish to find an interval [a, b] such that F(t) = on this interval. Suppose that 
F : [a, b] — > M is continuous and that u £ R such that: F(a) < u < F(b). Then, by the 
intermediate value theorem, 3 c £ [a, b] such that F(c) — u. Clearly, we wish to find a 
c such that F(c) = 0, that is, u = 0. This means that u — satisfies: 

F{a) = — + | < < F(b) = i + J (7.10) 



This implies that: 



3 v ' 6 t 3 



alb , 

o<0<7T+o ( 7 - U ) 



3 3 



Since in an expanding universe, 6% > 0, we need to choose [a, b] £ M. such that Eq. 
(7.11) is satisfied. Any choice for which a < 0, and b > will satisfy this relationship. 
In this interval: 

F ^ = J.+l= Q ^\ c \ = J. ( 7 - 12 ) 

Therefore, Fit) = does indeed equal zero at c £ [a, b] such that a < 0,b > 
and c — j-. Which means that at some time t = c, the expansion scalar was infinite, 
which indicates that there was a singularity at t = c. Note that, in the above analysis, 
t = is some arbitrary time, and not the time of the singularity. It can be thought 
of as the initial time, and 9i is the value of the expansion scalar at this time. The 
arguments given above revolve around the idea that the equation for 9(t) is integrated 
backwards upon which noting that there exists a time t S i ng < which corresponds to 
the singularity. It then so happens that the singularity occurs at t = t s < 0. 
A more formal version of the arguments given above could be given as follows. We 
essentially follow the methods in [33]. 

Lemma 1. Let F(t) = gjjj = j- + | be continuous on [a,b]. If F(a) < < F(b) V 
F(b) < < F(a) =£• 3 d, such that a < d < b for which F(d) = 0. 

Proof. Suppose that F(a) < < F(b). Since F(a) < 0, by the continuity of F(t), 
3 h £R such that F < on [a, h). Define: 

d = sup{/i : F{t) < on [a, h)} (7.13) 
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It is clear that d < b. One cannot have F(d) > 0, since this would mean that F > 
on some interval to the left of d, which is a contradiction based on our definition of 
d above. Additionally, we cannot have that d — b, so, d < b. We also cannot have 
F(d) < 0, since this would imply that there would be an interval [a, t), with t > d 
which would in turn imply that F(t) < 0. This, once again, contradicts the definition 
of d. Therefore, it must be that F(d) — 0. □ 

All of the aforementioned arguments point to a general theorem on singularities first 
given by Penrose and Hawking [18] [15] [17]. More or less, we can say that if matter 
as described by the energy-momentum tensor satisfies the strong energy condition, 
and there exists a a € M. + such that 8(t) > a, everywhere in the past of a specific 
spatial slice, then, there exists a singularity point where all geodesies end. As noted by 
Hervik and Gr0n, spacetimes can have singularities even if the strong energy condition 
is violated. If, however, we assume that the strong energy condition must be satisfied, 
then, for a viscous universe, Eq. (7.3) implies conditions on C for the existence of a 
past singularity: 



--K(ji + 3p) < -^kC 

^^(M + 3p) > C (7.14) 

Since this relationship was derived independent of any specific Bianchi type, one 
can conclude that for a past singularity to exist in an universe described by an energy- 
momentum tensor of the form Eq. (3.16), with an equation of state of the form Eq.(7.2), 
Eq. (7.14) must strictly hold. A diagram of an expanding congruence of geodesies from 
an initial singularity point is shown in Fig. (7.1) below. 
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Figure 34: E xpanding Congruence of Geodesies From an Initial Sin gularity Point 
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8 Conclusions 

In this thesis, we were interested in formulating a viscous model of the universe based 
on The Bianchi Type IV algebra. We first began by considering a congruence of fluid 
lines in spacetime, upon which, analyzing their propagation behaviour, we rigorously 
derived the famous Raychaudhuri equation, but, in the context of viscous fluids. We 
then went through in great detail the topological and algebraic structure of a Bianchi 
Type IV algebra, by which we derived the corresponding structure and constraint equa- 
tions. From this, we looked at the Einstein field equations in the context of orthonormal 
frames, and derived the resulting dynamical equations: The Raychaudhuri Equation, 
generalized Friedmann equation, shear propagation equations, and a set of non-trivial 
constraint equations. We showed that for cases in which the bulk viscous pressure is 
significantly larger than the shear viscosity, this cosmological model isotropizes asymp- 
totically to the present-day universe. We finally concluded by discussing The Penrose- 
Hawking singularity theorem, and showed that the viscous universe under consideration 
emerged from a past singularity point. 
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